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Summary. — Using a generalized parton model approach including spin and in- 
trinsic parton motion effects, and assuming the validity of factorization for large pr 
jet production in hadronic collisions, we study the azimuthal distribution around 
the jet axis of leading pions, produced in the jet fragmentation process. We identify 
the observable leading- twist azimuthal asymmetries, which are generated by all the 
physically allowed combinations of transverse momentum dependent (TMD) parton 
distribution and fragmentation functions. In particular, we show how one can isolate 
the Collins and Sivers contributions, and suggest a test of the process dependence 
of the Sivers function by considering the effect of color-gauge invariant initial and 
final state interactions. 

PACS 13.85.Ni - Inclusive production with identified hadrons. 
PACS 13.88.+e - Polarization in interactions and scattering. 



1. Introduction 



Transverse single-spin and azimuthal asymmetries in high-energy hadronic reactions 
have raised a lot of interest in the last years (see e.g. Ref. [1] and references therein). 
In particular, the huge spin asymmetries measured in the inclusive forward production 
of pions in high-energy proton-proton collisions, at moderately large transverse momen- 
tum, cannot be explained in the realm of leading- twist (LT) perturbative QCD (pQCD), 
based on the usual collinear factorization theorems. Out of the theoretical approaches 
proposed in order to account for these measurements, in the following we will adopt the 
so-called transverse momentum dependent (TMD) formalism, which takes into account 
spin and intrinsic parton motion effects assuming a pQCD factorization scheme. In this 
framework, single-spin and azimuthal asymmetries are generated by TMD polarized par- 
tonic distribution and fragmentation functions, among which the most relevant from a 
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phcnomcnological point of view are the Sivers distribution [2] and, for transversely polar- 
ized quarks, the Boer-Mulders distribution [3] and the Collins fragmentation function [4] 
(similar functions can be defined for linearly polarized gluons, see e.g. Ref. [5]). 

Azimuthal asymmetries in the distribution of pions inside a large transverse momen- 
tum jet are quite interesting obscrvables [6] and are presently under active investigation 
at the Relativistic Heavy Ion Collider (RHIC). In contrast to inclusive pion production, 
where several underlying competing mechanisms (mainly the Sivers and Collins ones) 
cannot be separated, and in close analogy with the semi-inclusive deep inelastic scat- 
tering (SIDIS) case, one could discriminate among different effects by taking suitable 
moments of these asymmetries. In principle, quark and gluon originating jets can also be 
distinguished, at least in some kinematic regimes. Furthermore, information on the size 
of TMD distributions and fragmentation functions can be gained, in a kinematic region 
in which they are still poorly known. This would be very helpful in clarifying the role 
played by the quark(gluon) Sivers distribution and by the Collins (-like) fragmentation 
function in the sizable single spin asymmetries observed at RHIC for single inclusive pion 
production. A similar analysis where transverse partonic motion was considered only in 
the fragmentation process, aimed at a study of the universality of the Collins function 
for quarks, was presented in Rcf. [7]. Our approach, named generalized parton model 
(GPM), has in principle a richer structure in the observable azimuthal asymmetries, be- 
cause intrinsic motion is taken into account in the initial hadrons as well. Preliminary 
data from the STAR collaboration at RHIC seem to support our model, since they report 
on a Sivers asymmetry for neutral pions larger than zero [8, 9] and compatible with our 
predictions. However, since factorization has not been proven in this case, but is rather 
taken as a reasonable phenomenological assumption, the validity of the scheme and the 
universality of the TMD distributions involved require a severe scrutiny by further com- 
parison with experiments. 

2. Theoretical framework 

Wc consider the process p^p — > jet+7r+ X, with one of the protons in a transverse spin 
state described by the four-vector S. We work in the pp cm. frame, where the polarized 
proton moves along the +Z cm direction, and define (XZ) cm as the production plane 
containing the colliding beams and the observed jet, with (f>j)x cm > 0. In this frame 
S = (0, cos 05, sin 0s, 0) and pj = pj T(cosh?7j, 1, 0, sinhr/j), where ryj = — log[tan(#j/2)] is 
the jet (pseudo)rapidity. We denote by z and k± n respectively the light-cone momentum 
fraction and the transverse momentum of the observed pion inside the jet w.r.t. the jet 
(fragmenting parton) direction of motion. Calculations have been performed within the 
GPM framework at leading-order (LO) in pQCD utilizing the helicity formalism. More 
details can be found in Ref. [6]. 

The final expression for the single-transverse polarized cross section has the following 
general structure: 

2dcr(0s, 4>n ) ~ der + dAcr sin 4> s + dai cos 4>^ + da 2 cos 24>" + dAtrf sin(0 s - 4>k ) 
(1) + dA<7+ sin(0 s +<f%) + dA(T 2 - sin(0 s ~2<ffi) + dAa+ sin(0 s + 2<ff), 

where (f>^ is the azimuthal angle of the pion three-momentum around the jet axis, as mea- 
sured in the fragmenting parton helicity frame [6]. In close analogy with the SIDIS case, 
in order to single out the different contributions of interest, we can define appropriate 
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where W(<frs , 4>k) 1S some appropriate circular function of <ps and . 

3. Phenomenology 

In this section we present and discuss some phenomenological implications of our 
approach in kinematic configurations accessible at RHIC by the STAR and PHENIX 
experiments. We consider both central (rjj = 0) and forward (t/j = 3.3) (pscudo)rapidity 
configurations at a cm. energy s/s = 200 GcV (different cm. energies, namely s/s = 62.4 
and 500 GcV, are also studied in [6]). 

For numerical calculations all TMD distribution and fragmentation functions are 
taken in the simplified form where the functional dependence on the parton light-cone 
momentum fraction and on transverse motion are completely factorized, assuming a 
Gaussian-like flavour- independent shape for the transverse momentum component. Con- 
cerning the parameterizations of the transversity and quark Sivers distributions, and of 
the Collins functions, we consider two sets: SIDIS 1 [10, 11] and SIDIS 2 [12, 13]. Further- 
more, for the usual collinear distributions, we adopt the LO unpolarized set GRV98 [14]. 
For fragmentation functions, we adopt two well-known LO sets among those available in 
the literature, the set by Kretzer [15] and the DSS one [16]. Our choice is dictated by 
the subsequent use of the two available paramctrization sets for the Sivers and Collins 
functions in our scheme, that have been extracted in the past years by adopting these 
sets of fragmentation functions. Since the range of the jet transverse momentum (the 
hard scale) covered is significant, we take into account proper evolution with scale. On 
the other hand, the transverse momentum component of all TMD functions is kept fixed 
with no evolution with scale. In all cases considered, kj_^ is integrated over and, since we 
are interested in leading particles inside the jet, we present results obtained integrating 
the light-cone momentum fraction of the observed hadron, z, in the range z > 0.3. 

We have considered first, for tt + production only, a scenario in which the effects of 
all TMD functions are over-maximized. By this we mean that all TMD functions are 
maximized in size by imposing natural positivity bounds. The transversity distribution 
has been fixed at the initial scale by saturating the Soffcr bound and then we let it evolve. 
Moreover, the relative signs of all active partonic contributions arc chosen so that they 
sum up additively. In this way we set an upper bound on the absolute value of any of the 
effects playing a potential role in the azimuthal asymmetries. Therefore, all effects that 
are negligible or even marginal in this scenario may be directly discarded in subsequent 
refined phenomenological analyses. See Ref. [6] for a more detailed discussion. 

As a second step in our study we consider, for both neutral and charged pions, only 
the surviving effects, involving TMD functions for which parameterizations are available 
from independent fits to other spin and azimuthal asymmetries data in SIDIS, Drell-Yan 
(DY), and e + e~ processes. In Fig. 1 we present, in the forward rapidity region, the quark 

generated asymmetry A s ^^ s ^™ ' , which comes mainly from the convolution between the 
transversity distribution and the Collins fragmentation function. The plots are obtained 
adopting the parameterizations SIDIS 1 and SIDIS 2. The Collins asymmetry for neutral 
pions results to be almost vanishing, in agreement with preliminary RHIC data [8, 9]. 
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Fig. 1. - The estimated quark Collins asymmetry for the p'p — > jet + ty + X process, obtained 
adopting the parameterizations SIDIS 1 and SIDIS 2 respectively, at yfs = 200 GeV in the 
forward rapidity region. The dotted black vertical line delimits the region xf ~ 0.3. 



The dotted black vertical line delimits the region xf ~ 0.3, with xf = %P}l/\/s, beyond 
which the SIDIS parameterizations for transvcrsity are extrapolated outside the Bjorkcn 
x region covered by SIDIS data and are therefore plagued by large uncertainties. For 
this reason, the two parameterizations for charged pions give comparable results only 
below this limit (notice the difference in scale between the two panels). A measurement 
of such asymmetries would be then very important and helpful in clarifying the large x 
behaviour of the quark transversity distribution. Furthermore, it turns out that in the 
central rapidity region the quark Collins asymmetries are practically negligible. 

In Fig. 2 we show, for both neutral and charged pions, the quark and gluon contri- 
butions to the Sivers asymmetry A s ™^ s , which cannot be disentangled, in the forward 
rapidity region as a function of PjT- The quark term is obtained adopting the SIDIS 1 
and SIDIS 2 parameterizations. The almost unknown gluon Sivers function is tentatively 
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The estimated quark and gluon contributions to the Sivers asymmetry for the p^p — > 
X process, obtained adopting the parametrization sets SIDIS 1 and SIDIS 2, at forward 
and = 200 GeV. The dotted black vertical line delimits the region xf ~ 0.3. 
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Fig. 3. - The Sivers asymmetry A B ^^ S for the p^p — > jet + ir + X process as a function of pjT, 
at fixed jet rapidity rj- s = 3.3, for the RHIC energy -Js = 500 GeV. 



taken positive and saturated to an updated version of the bound obtained in Ref. [17] by 
considering PHENIX data for the ir° transverse single spin asymmetry at mid-rapidity. 
Similarly to the case of the Collins asymmetry, the two parameterizations give compa- 
rable results only in the pj t region where the behaviour of the quark Sivers distribution 
is reasonable well constrained by SIDIS data (see again the dotted black vertical line). 

4. Test of the process dependence of the Sivers function 

In the GPM approach one applies the TMD distribution and fragmentation functions 
extracted from SIDIS, considering them to be universal. Very recently [18] the azimuthal 
distribution of leading pions inside jets has been studied allowing for process dependence 
of the (quark) Sivers function. This is referred to as the colour gauge invariant (CGI) 
GPM [19]. Namely, we have taken into account the effects of initial (ISIs) and final state 
interactions (FSIs) between the active parton and spectator remnants in the different 
scattering sub-processes. Since the details of such interactions depend on the particular 
partonic reaction considered, they render the Sivers function non-universal (see [19] and 
references therein). The oft-discussed case is the difference between the FSIs in SIDIS 
and the ISIs in DY scattering which leads to the prediction of an opposite relative sign 
of the Sivers functions in the two processes. This is considered to be a crucial test of 
our understanding of single spin asymmetries in QCD and still has to be confirmed by 
experiments. 

When applying similar reasoning to hadron production in pp collisions, typically the 
Sivers function has a more complicated colour factor structure since both ISIs and FSIs 
contribute. However, in the forward rapidity region the process under study is dominated 
by only one channel, qg — > qg, with the final quark identified with the observed jet. As 
a consequence, one finds that the predictions for the Sivers asymmetries obtained with 
and without inclusion of color gauge factors are comparable in size but with opposite 
signs [18], similarly to the DY case. Therefore the experimental observation of a sizable 
asymmetry could easily discriminate among the two approaches and test the universality 
properties of the Sivers function in hadronic reactions. Our results are shown in Fig. 3, 
where we plot A s ^ ^ s as a function of the jet transverse momentum p 3 T at r/j = 3.3, for 
the RHIC energy y/s = 500 GeV, integrated over fej_ w and z (z > 0.3). The predictions 
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labeled SIDIS 1 and SIDIS 2 are similar in the intermediate p^r < 5-5 GeV region 
(corresponding to xp < 0.3). This region is then optimal to test directly the process 
dependence of the Sivers function. Note that for y/s = 200 GeV the behavior of our 
estimates would be similar to that shown in Fig. 3, gaining almost a factor of 2 in size. 
However the range of p 3 T covered would be narrower (py? < 6.5 GeV) and with xf < 0.3 
now corresponding to p^r < 2.5 GeV. 

Finally, as a natural extension of this analysis, one can consider single inclusive jet 
asymmetry in pp scattering, which is described solely by the Sivers function. The results 
obtained for A s ^ ^ s (not shown) look almost indistinguishable from the case of neutral 
pion-jct production (central panel of Fig. 3). 

5. Conclusions 

Wc have presented a study of the azimuthal asymmetries measurable in the distribu- 
tion of leading pions inside a large-pr jet produced in single-transverse polarized proton 
proton collisions for kinematic configurations accessible at RHIC. To this end, we have 
adopted a generalized TMD parton model approach with inclusion of spin and intrinsic 
parton motion effects both in the distribution and in the fragmentation sectors. In con- 
trast to inclusive pion production, where the Sivers and Collins mechanisms cannot be 
separated, the leading-twist azimuthal asymmetries discussed above would allow one to 
single out the different effects by taking suitable moments of the asymmetries. This will 
give us the opportunity of testing the factorization and universality assumptions, and of 
gaining information on the size and sign of the various TMD functions in a kinematic 
region not covered by SIDIS data. 
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